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ABSTRACT 

We carry out fully 3-dimensional simulations of evolution from self-similar, spherically sym- 
metric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As 
a result of the radial orbit instability, the haloes which grow from such initial conditions are 
triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately 
self-similarly in time. In all cases they have power-law density profiles and near-constant ve- 
locity anisotropy in their inner regions. Both the power-law index and the value of the veloc- 
ity anisotropy depend on the similarity index of the initial conditions, the former as expected 
from simple scaling arguments. Halo structure is thus not "universal" but remembers the ini- 
tial conditions. On larger scales the density and anisotropy profiles show two characteristic 
scales, corresponding to particles at first pericentre and at first apocentre after infall. They are 
well approximated by the NFW model only for one value of the similarity index. In contrast, 
at all radii within the outer caustic the pseudo phase-space density can be fit by a single power 
law with an index which depends only very weakly on the similarity index of the initial con- 
ditions. This behaviour is very similar to that found for haloes formed from ACDM initial 
conditions and so can be considered approximately universal. 
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1 INTRODUCTION 

Galaxies are believed to form from gas condensing at the centres 
of massive dark matter haloes as these grow by collapse and ag- 
gregation from weak density fl uctuations emerging from the early 
Universe ( jWhite & Rees||l978i) . 

The earliest theoretical insights into the formation and evo- 
lution of dark matter haloes w ere provided by the spherical infall 
model of iGunn & Got^ h972h and lGotJ ( Il975l) . In this model, an 
isolated overdensity in an otherwise unperturbed Einstein-de Sitter 
universe first expands with the Hubble flow, then turns around and 
collapses. Surrounding material continues to fall onto the object un- 
til its mass greatly exceeds that of the originally perturbed region. 
Asymptotically, a power-law density profile is established with 
p oc r~^'^^ . The late- time structure of this model i s a similarity 
solutio n whose structure was worked out in detail bv lBertschingen 
( 119851) . More general similarity solutions where the initial mass 
perturbation scales with enclosed mass as SM/M oc M"" were 
presented by iFillmore & GoldreichI j 19841) who showed that these 
produce haloes with p oc r~"', where 7 = 2 for < e ^ 2/3 
and 7 = 9e/(l + 3e) for e ^ 2/3. The change in behaviour at 
e = 2/3 reflects the fact that all orbits are purely radial in these 
models. No superposition of radial orbits can self-consistently pro- 
duce a power-law shallower than r~^ at small r. Further gener- 
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alisations to include randomly oriented rosette orbits with a scale- 
free eccentricity distribution resulted in spherically symmetric sim- 
ilarity solutions for whic h 7 = 9£/(l + 3e) for all e > 
(White &Zaritsk v 1992; Sikivie et alj|l995l . Il997l : [Nussedl200ll : 
iMacMillan et alj|2006l) . In such models the inner density structure 
of the final haloes reflects the scaling properties of the initial con- 
ditions. 

In contrast, N-body simulations have shown for more than a 
decade that dark matter haloes formed from fully 3-dimensional, 
cosmologically realistic initial conditions do not have pure power- 
law density profiles. Rather, the logarithmic slope of simulated pro- 
files changes slowly but continuously with radius. In addition, the 
shape of these profiles is almost independent of halo mass, of cos- 
mologic al parameters, and o f the power spectrum of initial fluc- 
tuations dNavarro et al.lll997D . The most popular representation of 
this "universal" shape is the NFW model wh ich behaves as 1/r i n 
the inner regions and as l/r'* in the outskirts dNavarro et al J 19971) . 
Even haloes formed by monolithic collapse in the first nonlinear 
phases of hot and warm dark matter models (HDM and WDM) are 
well represented by NFW fits, showing that hie rarchical growth i s 
not required to p roduce this universality (see e.g. lHuss et alj ( ll999f) : 
IWang & White! ( 120091) ). Recent N-body simulations with signifi- 
cantly increased numerical resolution have found small but signif- 
icant deviations from NFW shape which depend systematically on 
halo mass (Havashi & White 2008; Navarro et al. 2010). Thus halo 
density profiles are not truly universal. Nevertheless, these varia- 
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Figure 1. Projected density maps of the haloes formed for e = 0.4 (left) and e = 0.8 (right). In each case, the long axis of the bar is horizontal, its short 
axis is vertical, and the region plotted is a slice of thickness 2router caustic and side 4routor caustic where r^vLtav caustic is the radius of the outer caustic 
(first apocentre after the turnaround) in the spherical similarity solution. The thick and thin green circles have radii of router caustic and 2routcr caustic 
respectively. In units of the outer-caustic radius, the bar length is approximately "universal", although time-dependent features are still visible in these plots, 
for example, in the lack of left-right symmetry. 



tions are much smaller than those predicted by the similarity so- 
lutions, and the profiles can be considered universal to a good ap- 
proximation. 

The differences in behaviour between similarity solutions and 
numerical simulations reflect the fact that the former enforce spher- 
ical symmetry and a potential which varies smoothly with time ac- 
cording to the similarity scaling, while the latter involve strongly 
time-dependent and fully three-dimensional potential fluctuations 
which exchange energy and angular momentum between differ- 
ent parts of the system through the classical "violent relaxation" 
mechanism JLvnden-Belll 1 19670 . Even though extensions of the 
original similarity solutions allow non-radial orbits with a variety 
of eccentricity distributions jWhite & Zaritskvll 19921 : Isikivie et al.l 
1 1 99511 1 9971 : lNusserll200ll ; lMacMillan et al.ll2006l) . thev retain both 
spherical symmetry and strict similarity scaling, and so exclude 
violent relaxation and any possibility for it to drive convergence 
towards a universal (i.e. e-independent) nonlinear structure. Large 
and time-dependent potential fluctuations are required to produce 
such convergence. These occur naturally in hierarchical assembly 
models but also in other situations, for example, during the mono- 
lithic, quasi-ellipsoidal collap se of the first generation of haloes in 
HDM or WDM cosin ogonies l lZel'Dovichlll970l : ISheth et alj200ll : 
IWang & Whit3l2009l) . 

It has long been known that spherical equilibria dominated 
by radial orbits are violently unstable and evolve on a few dy- 
namical times i nto strongly eU i psoidal bars with significant non- 



radial motions jAnto noV 1973 



Polvachenko & Shukhman:.1981 



198J: 



Meza & Zamoran 



HenonI 1 19731: IPolvachenkd 1 198 ll: 
Barnesll 19851 : iMerritt & Aguilaj 



1997). Th e original sim i larity solutions 



of iFillmore & GoldreiclT l l 19841) and lSertschingeil ( 1 19851) are thus 
not viable models for the formation of real systems - the slightest 
non-spherical perturbation of their quasi-equilibrium inner regions 



causes rapid evolution to an entirely different nonlinear structure. 
A number of authors have noted situations where the radial or- 
bit instability gives rise to haloes with NFW- l ike density profiles 
JHenriksen & Widrowlll99^: IHuss et al. | l 19991 : ISarnes et alj|2005l : 
iMacMiUan et alJl2006inBellovarv et alj|2008l) . There are thus (at 
least) two possibilities for the long-term evolution of structure from 
(almost) spherically symmetric, self-similar linear initial condi- 
tions. Either it may approach a non-spherical self-similar solution, 
which would then have a potential which changes smoothly in time 
and no violent relaxation, or the strongly time-dependent behaviour 
may continue indefinitely, allowing violent relaxation to rearrange 
material in the inner regions. In the former case an NFW-like uni- 
versal profile is only consistent with similarity scaling for e = 1/6, 
whereas in the latter case an NFW-like profile could, in principle, 
be maintained at all times. 



In the present paper we investigate these issues by simulat- 
ing evolution from spherically symmetric, self-similar, linear ini- 
tial conditions for a var iety of values of e. As we showed in 
IVogelsberger et alj 1 120091) particle noise leads to the rapid onset of 
the radial orbit instability in such simulations, so that their later 
nonlinear evolution is strongly non-spherical. Here we show that 
while some chaotic time-dependence remains at late times, evo- 
lution is nevertheless approximately self-similar. The inner struc- 
ture of the haloes remembers the initial conditions from which they 
formed, depending on e in the same way as in the spherical similar- 
ity solutions. In Section[2]we present our simulations and study the 
formation of ellipsoidal "bars". In Section [3] we study how these 
bars affect the velocity anisotropy and density profiles in the inner 
halo. The transition between the infall and quasi-equilibrium re- 
gions, and the inner and outer scales which define it are discussed 
in Section|4l Finally in Section[5]we demonstrate that, in contrast to 
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Figure 2. Left panel: Halo axis ratios as a function of radius for a variety of similarity scaling parameters. Different colours correspond to different values 
of the scaling parameter e. For each r, axis ratios are estimated from the principal values of an inertia tensor calculated for all particles within an ellipsoid of 
mean radius r. Dashed lines denote minor-to-major axis ratios while solid lines denote intermediate-to-major ratios. For each simulation, radii ai'e normalised 
by that of the outer caustic in the corresponding similarity solution. The left panel shows results for the final time a = 1. In the right panel these are replotted 
as thick lines and compared with results for a = 0.5 (the thin lines). 



the velocity anisotropy and density profiles, the profile of pseudo- 
phase-space density is almost universal in these models. 



2 RADIAL-ORBIT INSTABILITY: BAR FORMATION 

The similarity solutions assume an Einstein-de Sitter universe in 
which the linear mass perturbation SAIi within a sphere containing 
unperturbed mass Mi, initially satisfies 



finding the inner profile to vary with the initial scaling properties: 



SMj 



1.0624 



Ma 



(1) 



where e is a scaling index and Mq is a reference mass taken to 
be the mass within the turnaround radius (rta) at the initial time. 
The parameter e is restricted to positive values in order to ensure 
that more distant mass shells turn around and fall back later than 
inner ones. The mass within the shell that is just turning around 
and the physical radius of this shell at turnaround then scale with 
time as Mta oc t^^^" and rta oc ^2/3+2/9^ ^^ jj^^j jj^^ mewa mass 
density contained within the turnaround radius is always 5.5 times 
the cosmic mean (the critical density) and so satisfies pta oc t^^ . 

Since these initial conditions are spherically symmetric, the 
system remains spherically symmetric at later times (if instabilities 
are ignored/suppressed) and the particles on each spherical shell 
all execute identical though differently oriented radial orbits. The 
scale-free nature of the initial conditions ensures that the orbits ex- 
ecuted by different shells are i dentical when scaled to their individ- 
ual turnaround radii and times. iGurm & GottHl972h showed that at 
radii much smaller than rta this results in a stable halo with den- 
sity profile p oc r~^'^ when e = 1, the only case they considered. 
iFillmore & GoldreichI ( Il984f ) generalised this to other values of e. 



p(r) oc 



-9e/(1+36) 






(2) 



The change of behaviour at e = 2/3 is due to the fact that no self- 
consistent spherical equilibrium system built from radial orbits can 
have a power-law slope shallower than —2, because every parti- 
cle is constrained to pass through the centre of the system once 
per orbit. The full structur e of the simil a rity so lution for the case 
e = 1 was worked out by iBertschingen l ll985l) . who showed that 
the power-law inner structure breaks in the transition to the infall 
regime, and that the confinement of non-zero phase-space density 
to a 3-dimensional sheet results in a series of sharp spherical caus- 
tics which are superposed on the inner power-law density profile 
and are located at any given time at the positions of particles cur- 
rently passing through apocentre. 

In this paper we simulate evolution from initial condi- 
tions which obey the sph erical similarity solutions, exactly as in 
IVogelsberger et alj 1 120091) . but for a variety of values of e. We use 
a fully three-d imensional N-body solver, a version of Gadget-2 
l lSDringe]||2005l) . and we allow particle noise to drive the radial or- 
bit instability. This happens in the first few expansion factors so 
that the later stages of the simulations all contain fully developed 
ellipsoidal "bars". All our simulations follow evolution over a fac- 
tor of 1000 in time corresponding to expansion of the background 
cosmology from a — 0.01 and until a = 1. The softening length 
is kept fixed in comoving coordinates at 0.00025rta(a = 1). The 
final halo is represented by about 256'' /2 ~ 8.4 x 10® particles 
within rta. 

In Fig. [T] we show images of the bar in the final state of two 
of our simulations, e — 0.4 on the left and e = 0.8 on the right. In 
each case the bars extend slightly beyond the position of the outer 
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Figures. Left panel: Velocity anisotropy /3(r) = 1 — erf /a'^ as a function of radius r, normalised by the turnaround radius rta(i)- Here ct^C'') and crf{r) sue 
the 1-dimensional radial and transverse velocity dispersions averaged over spherical shells. Different colours denote simulations with different values of e as 
indicated. Solid lines aie for a = 1 and dashed hnes for a = 0.5. Beyond the outer caustic, all particles are falling in for the first time, but they nevertheless 
have nonzero radial and tangential velocity dispersions because the quadrupole moment of the ellipsoidal inner regions induces variations in infall velocity 
and infall direction at each r. Both dispersions jump dramatically as the outer caustic is crossed and the anisotropy is slowly varying and relatively small in 
the inner quasi-equilibrium regions. The softening radius is too small to affect the velocity dispersions over the radial range plotted. Right panel; Mass density 
averaged in spherical shells as a function of radius, again normahsed by rta. Line colours and types coiTespond to those in the left panel. The density is given 
in units of the current critical density and is multiplied by (r/rta)^ to suppress the dominant dependence and to allow easier comparison of the curves. In 
both panels the structure is very similar at the two times considered. The dependence on e is very clear in the density profiles and appears significant also in 
the anisotropy profiles. 



caustic of the corresponding spherical similarity solution, which we 
indicate with a green circle. The location of the actual outer caustic 
corresponds to the transition in colour from blue to red. In all our 
simulations we find bars with semi-major axes (delineated by the 
outer caustics) which are about L2 times the radius of the spherical 
caustic. Thus measured in units of the latter, our bar lengths appear 
universal. In both the cases shown there are clear left-right asym- 
metries in the inner structure of the bars (the yellow regions). Such 
asymmetry is a clear indicator of time-dependent behaviour - any 
exact non-spherical similarity solution is expected to be left-right 
symmetric. 

We show the axis-ratios of equidensity surfaces of our final 
bars as a function of semi-major axis in the left panel of Fig. [2] 
These ratios were obtained by calculating the moment of inertia of 
all particles within an ellipsoidal surface, deriving the eigenvalues 
and principal axes, and then iterating until the shape and orien- 
tation of the bounding ellipsoid are consistent with the orientation 
and relative axis lengths inferred from the moment of inertia. These 
ratios are thus cumulative, referring to all material interior to the 
quoted radius rather to an ellipsoidal shell at this radius. The axis 
ratios take their most extreme values just inside the outer caustic, 
with all simulations giving minor-to-major values close to 0.2 and 
intermediate-to-major values close to 0.35. They approach unity 
rapidly at larger radii, and also rise steadily though more slowly 
towards smaller radii. There is no clear trend as a function of e, 
suggesting that the exact values may be stochastically determined 
and perhaps also time-dependent. We test this in the right panel by 
replotting the curves using a thick line-style and comparing them 
with results at a = 0.5, represented by thin lines. The overall pat- 



tern is very similar at the two times. There is some indication that 
the deviation of individual models from the mean is coherent be- 
tween the two times, but there are also substantial variations. Time- 
dependent effects appear to be influencing the measurements sig- 
nificantly, and these plausibly have a relatively long time coher- 
ence. In consequence it is not possible to decide whether the appar- 
ent trends with e are real, or just reflect unrelated time-dependent 
variations. The overall similarity of all the curves at both times sug- 
gest that the shape behaviour is approximately both self-similar and 
"universal". 



3 DENSITY AND VELOCITY ANISOTROPY PROFILES 

The central bar torques infalling particles, inducing non-radial mo- 
tions which transform the purely radial orbits of the original simi- 
larity solution into fully three-dimensional orbits. This alters both 
the density and the velocity dispersion structure of the nonlinear 
regions, as we now illustrate. A velocity anisotropy parameter is 
conventionally defined as 



/3(r) 



2 / 2 



(3) 



where ar{r) and crt(r) are 1-dimensional radial and transverse ve- 
locity dispersions averaged over spherical shells of radius r. The 
left panel of Fig. [3] shows /?(r) at two different times, a = 0.5 
(dashed lines) and a — 1 (solid lines). The right panel shows 
{r/rti,)'^p(r)/pb at the same two times. The two quantities are 
evaluated using the same set of logarithmically spaced spherical 
shells. Neither profile changes systematically with time when plot- 
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ted against r/rta, consistent witli the expectations of self-similar 
evolution. 

Tlie splierically-averaged density profiles in the right panel of 
Fig. [3] show three distinct regimes. At radii larger than the outer 
caustic, only one (infalling) matter stream is present. Here, there 
is essentially perfect agreement between a = 0.5 and a = 1 in 
all cases, and the effective power-law index of the profile is signifi- 
cantly greater than —2 and similar for all e. Only the position of the 
outer caustic (in units of rta) varies significantly with e. Between 
this outer caustic and the typical radii of particles passing their first 
pericentre after tumarounqj, the density profile has a mean effec- 
tive index which is similar for all e and substantially less than —2. 
The behaviour in this regime is quite irregular, however, and there 
are substantial differences between a = 0.5 and a = 1 in several 
ca ses. Examination of phase -space plots (see, for example, Fig. 8 
of IVogelsberger et alj j2009l) ) suggests that this time-dependence 
is driven by large-scale irregularities which grow around the first 
apocentre after turnaround but are washed out by phase-mixing at 
later times. Inside the radius of first pericentre, the profiles become 
more regular and agree well between the two times shown. To a 
good approximation, they are power laws that are well fitted (as 
r — >• 0) by 



p{r) OH r ' ^ , e > 0. 



(4) 



This is the value expected for a pe rfect similarity solution 
in the presence o f nonr a dial motions JWhite & ZaritskyI 1 19921 ; 
ISikivie et alj 1 1991 1 19971 : iNusseJ l200ll : iMacMillan et alj 1200^ . 
demonstrating that the time-dependent behaviour seen at interme- 
diate radii is not sufficient to destroy the similarity scaling and 
enforce "universal" structure. Note that both the steep slope over 
the decade immediately inside the outer caustic radius (which is 
roughly at the conventional virial radius) and the break to a shal- 
lower slope at smaller radii are qualitatively similar to the "uni- 
versal" behaviour encapsulated by the NFW profile. However, the 
inner slope matches the NFW value only for e — 1/6, smaller than 
any of the values tested here. 

The left panel of Fig.[3]shows that velocity anisotropy profiles 
are similar for all values of e, exhibiting distinct behaviour in the 
same three regimes seen in the density profiles. At large radii, par- 
ticles are falling in for the first time on very nearly radial orbits, 
gradually gaining angular momentum as they feel the quadrupole 
moment of the central bar. As they cross the outer caustic they 
are mixed with particles which have already passed pericentre, and 
/3(r) drops to much smaller values. Inside the radius of first peri- 
centre, /3(r) is almost constant in the quasi-equilibrium region, de- 
clining slightly towards the centre. In this region there appears to 
be a weak but significant dependence of anisotropy on similarity 
index, with larger values of e giving rise to slightly more radially 
biased velocity dispersions. This behaviour is a consequence of the 
dependence of the inner density profile on e. For larger e the in- 
ner mass distribution is more strongly centrally concentrated, so 
that its quadrupole moment become less significant relative to the 
monopole and less angular momentum is induced in the orbits of 
infalling particles. 

The left panel of Fig. [3] suggests that in the inner part of the 
halo the (3 parameter may approach an asymptotic value which de- 
pends on the value of e. [Hansen & Moord 1 120061) have proposed 
that the equation /3(r) = -0.2(d(lnp)/d(lnr + 0.8)) relates ve- 
locity anisotropy and density profile slope in the inner parts of 



^ We investigate these scales in more detail below. 
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Figure 4. The radius of first apocentre after turnaround (i.e. the outer caus- 
tic radius) in units of the current turnaround radius as a function of e in the 
spherically symmetric similarity solutions with radial orbits. 
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Figure 5. Histograms of first pericentre distance, measured directly from 
our simulations, in units of the fiducial current outer caustic radius given in 
Fig. O In these units the distribution of first pericentre distance is almost 
independent of e, peaking at about 0.1. 



dark matter haloes (see also e.g. lNavarro et alJoOlOl) ). Using Eq.|4] 
which is indeed valid in the inner part of our haloes, the above re- 
lationship reduces to 



0.2 



9e 



0.8 , e > 0, 



which gives /3 



(1 + 36) 

{0.17,0.2,0.23,0.25,0.29} for e 



(5) 
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Figure 6. Using the radius of tlie outer caustic as a cliai'acteristic scale maJces the nonlinear stracture of haloes appear almost independent of e. Left and right 
panels show the effect of using this scaHng on the velocity anisotropy and density profiles, respectively. In the right panel, the j/-axis is furthermore multiplied 
by an e-dependent power law of radius chosen so that the inner profile appears flat and equal to unity (see Eq.|4j- Solid lines show the a = 1 profiles, whereas 
dashed lines show the profiles at time a = 0.5. 



{0.4, 0.5, 0.6, 0.7, 0.9}. These values are in good agreement with 
those infeffed from the left panel of Fig. [5] even though our haloes 
form in a quite different (and "nonuniversal") manner from CDM 
haloes. 



4 THE INNER AND OUTER SCALES : FIRST 
PERICENTRE AND OUTER CAUSTIC RADII 

In the previous section, we found approximate similarity behaviour 
for nonspherical haloes formed by infall from power-law and al- 
most spherically symmetric initial con ditions. As in the full y spher- 
ically symmetric similarity solution of lBertschingeg jl985l) . the re- 
sulting nonlinear density profile is a power law only at sufficiently 
small radii, but in our case the transition between the infall and 
quasi-equilibrium regimes is more complex than when the orbits 
are purely radial. As we have seen, there appear to be two distinct 
scales: an iimer scale at around O.Olrta and an outer scale at around 
O.lrta. In the last section we asserted that the former can be associ- 
ated with the first pericentric passage of infalling particles, and the 
latter with the following apocentric passage, which occurs at the 
outer caustic. Here we support these assertions by a more careful 
study of the radial position of these orbital turning-points. 

In the spherical similarity solution, the outer caustic is the bor- 
der between the one- and three-stream regions. Outside this radius 
all particles are falling in for the first time, while within it many 
of them have passed through the centre at least once. In Fig.|4]we 
show how the radius of this outer caustic, evaluated directly from 
the similarity solution, increases with e. Comparison with the right 
panel of Fig. [3] shows that the sharp change in profile shape which 
delineates the inner boundary of the infall regime in our simula- 
tions is indeed very close to this radius and depends on e exactly as 
predicted. 

In Fig. [3] all radii are expressed in units of the turnaround ra- 



dius. Although this is the natural scale for models with spherically 
symmetric initial perturbations of the kind studied here, it turns out 
that scaling to the outer caustic radius (i.e. the position of the first 
apocentre after turnaround) results in greater uniformity of the in- 
ner nonlinear structure of our haloes as a function of e. We show 
this in Fig.[6]which may be compared with Fig. [3] Note that in the 
right panel of Fig.|6l the outer transition occurs slightly beyond the 
nominal outer caustic radius. This is because we have rescaled us- 
ing the value for the spherically symmetric similarity solution, as 
plotted in Fig. [4] This is slightly smaller than caustic radius at the 
end of our numerical bars, as is clearly visible in Fig.[T] 

Next, we consider the inner scale at around O.Olrta (see 
Fig. [3}. We follow the trajactories of all simulation particles and 
record the radial position of their first pericentric passage. For all 
particles which pass first pericentre between a = 0.5 and a = 1, 
we follow their radial position r along their orbit and record its first 
minimum, as well as the time when this minimum occurs. The latter 
allows us to calculate the position of the outer caustic (according to 
the spherical similarity solution) at the time of pericentric passage, 
and thus to measure the ratio of these two lengths. A histogram 



01 fpcvi/'o' 



utcr caust 



ic is presented in Fig.|5] There is considerable 



scatter in this ratio because pericentric radius depends strongly on 
the angle between the infall direction of a particle and the long axis 
of the bar. It is very small for infall along one of the principal axes, 
and it maximises at angles well away from any of these axes. Never- 
theless, there is a well-defined "typical" pericentric distance which 
is ~ O.lrouter caustic for all valucs of e. This agrees well with the 
inner transition scale of the /3 and density profiles in Figs.[3]and|6l 



5 PSEUDO-PHASE-SPACE DENSITY PROFILES 

A property of dark matter haloes with rather intriguing charac- 
teristics is the pseudo-phase-space density. This is defined as, 
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Figure 7. Left panel: Pseudo-phase-space density profiles for various e evaluated at a = 0.5 (dashed) and a = 1 (solid). The radial coordinate is scaled with 
rta and the pseudo-phase-space density with /ta = Pb^^ /'''ta.- With these scalings, the profiles do not change significantly with time, and, remai'kably, ai'e 
very similar and are close to a power law within the outer caustic for all e, even though though they differ in the infall region. There is no obvious feature at the 
radius of first pericentre. Right panel: Pseudo-phase-space density profiles of the inner parts of our haloes after reseating the horizontal and vertical axes by 
'"ta/^'outor caustic and by (router caustic/rta)^ respectively. This scaling brings the overall profiles into quite close agreement within the outer caustic, as 
can be seen in the inset. The slope of the inner profile nevertheless depends weakly on e. In the main panel we have additionally scaled by r^**"'', the inverse 
of the behaviour expected theoretically for e = 1. Dashed lines show the pseudo-phase-space density profiles at a = 0.5. The thin straight lines in the lower 
left indicate the slopes of the analytic predictions of Eq.|6] 



p/(T'\ where p{r) and air) are the mass density and the (three- 
dimensional) velocity dispersion averaged over a shell of radius r. 
Even though neither p(r) nor oij-) is itself close to a power law, 
this particular combination is very close to a power law for haloes 
formed from ACDM initial conditions, and, remarkably, its power- 
law index is very sim ilar to that found in t he spherically symmetric 
similarity solu tion of Bertschingei 1 19851) . This curiosity was first 
pointed out bv lXavlor & Navarrd d200lh . and since that time it has 
been investigated by many authors who have found the behaviour 
to be quite robust, to extend over more than three decades in ra- 
dius, and to hold also for other kinds of initial co nditions (see, for 
example, iLudlow et aUbOld ; iNavarro et ai]|201(ll) ) . 

In the framework of a true similarity solution with nonradial 
motions the inner density profile must obey Eq.|4]for all e > and 
/3(r) must become constant at small radii. Since the circular ve- 
locity and the velocity dispersions must have the same scaling be- 
haviour, the pseudo-phase-space density is also a power law which 
is easily verified to be 



-3(2+3e)/(2(1+3£)) 



for all e 



(6) 



For th e original similarity solutions of Ipillmore & GoldreichI 
( 119841) which had purely radial orbits, this beh aviour holds only for 
e > 2/3. The specific case e = 1 studied bv lBertschingeii l ll985r) 
gives p/g^ cc r"-'^/*. 

We have directly evaluated the pseudo-phase-space density 
profiles of our simulated haloes using logarithmic bins in radius and 
subtracting the mean radial motion before evaluating the velocity 
dispersion. We show the results in Fig.|7] In the left panel the radial 
coordinate is normalised by the turnaround radius and the pseudo- 



phase-space density by the characteristic value /ta = ptt'^/rta- 
The very good agreement between the curves for a — 0.5 and for 
a = 1 demonstrates that the pseudo-phase-space density evolves 
self-similarly in time. This is no surprise, given that we have al- 
ready seen good scaling for the density and velocity anisotropy 
profiles. More surprising is the fact that inside the outer caustic 
radius the profiles show very little dependence on e, either in slope 
or in amplitude. Furthermore, they are all good approximations to 
a power law, and there is no obvious feature near the radius of first 
pericentre despite the very strong features seen at this radius in 
Fig.H] Once again scaling to the radius of the outer caustic gives an 
even better overlap since it matches the break at large radii, as can 
be seen in the inset in the right panel. 

We focus on the innermost part of our haloes (the regions in- 
side the first pericentre) in the main plot of the right panel of Fig. [7] 
We again use the outer caustic radius to scale the radii and pseudo- 
phase-space densities, and we additionally multiply the pseudo- 
phase-space density by r^^'^ to take out the dominant trend. The 
power-law slopes expected for exact similarity scaling are indicated 
by thin straight lines and are seen to be a good but not perfect fit to 
the simulated behaviour. The differences are small enough that they 
can probably be ascribed to residual time-dependent behaviour, an 
interpretation which is supported by the differences between the 
a = 0.5 and a = 1 curves for each e. 

Thus, unlike the velocity anisotropy and the density profiles, 
the pseudo-phase-space density profile is close to a power law over 
the full nonlinear extent of our haloes and it depends only very 
weakly on e. Thus, to a good approximation it can be considered 
"universal". 
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6 CONCLUSIONS 

We have shown that ahhou gh the classic spher i cally sym- 
metric similarity s olutions of iFillmore & GoldreichI ( 1 1984) and 
ISertschingeii ( 119851) are violently unstable to the radial orbit insta- 
bility, evolution from the initial conditions they presuppose gives 
rise to ellipsoidal, bar-like haloes which nevertheless grow in an 
approximately self-similar way. 

The nonlinear structure of these objects shows two character- 
istic radii. The outer caustic separates the infall and multistream re- 
gions and occurs at approximately the same position as in the orig- 
inal spherically symmetric similarity solutions. The second charac- 
eristic radius is about an order of magnitude smaller and occurs at 
the typical first pericentre distance of the infalling particle stream. 
It can be considered the outer edge of the quasi-equilibrium region. 
Relative to the turnaround radius, both radii increase with the value 
of the similarity parameter e. 

Both the density and velocity anisotropy profiles show strong 
features at the first pericentre radius. Within this radius the den- 



sity profile is approximately a power law, p oc 



with 7 



9e/(l + 3e), and the velocity anisotropy is approximately constant, 
/3 = 0.2(7 — 0.8). At larger radii the density profile becomes sub- 
stantially steeper and the velocity anisotropy rises steeply. 

Despite the strong features in the density and velocity 
anisotropy profiles at the first pericentre radius, the pseudo-phase- 
space density profiles of all our haloes are close to power laws 
all the way out to the outer caustic radius. Furthermore they de- 
pend only very weakly on e and are similar to the profiles of haloes 
formed from ACDM initial conditions. Thus this profile seems re- 
markably "universal". 

Our simulations were carried out at the Computing Centre 
of the Max-Planck-Society in Garching. RM thanks French ANR 
OTARIE for support. We thank Niayesh Afshordi, Ed Bertschinger, 
James Binney, Jacques Colin, Aaron Ludlow and Julio Navarro for 
useful discussions. 
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